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ABSTRACT 

The objective of this thesis was to develop a computer 
program for use in evaluating data acquired from the Depart- 
ment of Aeronautics Cessna 310H model aircraft. The program 
was designed for use with the IBM 360 computer, in conjunc- 
tion with a CALCOMP plotter. The program format is very 
similar to that of IBM's Scientific Subroutine Package 
routines; in addition, the program provides the student with 
a detailed set of instructions for loading the data. 

Using information taken from the engine curves for the 
aircraft, the program evaluates input data, determines the 
power required for flight at altitude, reduces the values 
to equivalent sea level values and outputs a graph of power 
versus velocity for level flight. Lift and drag coefficients 
are calculated as functions of aircraft weight, velocity, 
and power required, and are plotted on a second graph. 
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I. INTRODUCTION 



The laboratory course is one of the most demanding for 
the graduate student because of the extremely time consuming 
process of reducing raw data. Often the student may find 
himself neglecting important course work because of the 
large amount of work required to reduce laboratory data by 
hand to meaningful results. This is, of course, an unsat- 
isfactory situation, but is one which seems to be the rule 
rather than the exception for many students. 

The Department of Aeronautics Cessna 310H aircraft 
flying laboratory, which supports the Flight Evaluation 
Techniques courses, AE 3321 and AE 3322, at the Naval Post- 
graduate school, is typical of such time consuming labora- 
tories. With the incorporation of a photo-panel data 
acquisition system in the aircraft came the capability to 
collect large amounts of in-flight data which require large 
amounts of the student's time for evaluation. 

There appeared to be a need for some sort of procedure 
which would reduce the student workload in this laboratory 
by relieving him of the burden of hand reducing the large 
volume of raw data and hand plotting the results. A program 
which would require only that the student assemble the data 
in a proper sequence and input it to a computer, which in 
turn would evaluate the data and output tabular and graphical 
displays of the results, would be extremely beneficial in 
lightening the student workload. Such a program would free 



the student, enabling him to spend time evaluating his 
results instead of tabulating them. With this need in mind 
a computer program for the reduction of in-flight data for 
the Cessna 310H aircraft was written. 

CESSNA was designed as a student aid, with the intent 
that it be self-explanatory, containing all information and 
instructions necessary for the student £o run the program. 
Parts II and III of this thesis may be used as a user's 
manual in the event more detailed information is required, 
or if modification and/or additions are desired. The program 
is written in FORTRAN language but knowledge of the language 
is not required, other than the ability to type data cards 
and a job card in proper format; this format information is 
readily available at the computer center o The card deck is 
Complete. The student needs only to add his data cards to 
the end of the deck and it is ready to run. 
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II. PROGRAM ORGANIZATION AND OPERATION 



A. GENERAL ORGANIZATION 

As shown in Fig. 1, the complete program consists of 
two primary sections and two subroutines: 

1. Aircraft Power Requirements : 

2. Lift and Drag : 

3 . Interpolation Subroutine : 

4. Sorting Subroutine . 

When program execution begins, aircraft engine data are 
read in for internal use in the program. Following this, 
student data concerning aircraft weight, altitude, velocity, 
RPM and manifold pressure, and temperature, are read in for 
evaluation. 

CESSNA uses the engine data and the student's RPM and 
manifold pressure data to determine shaft horsepower required 
for level flight at the operating altitude for the input 
velocity range. Horsepower, along with velocity and weight, 
is then used to determine lift and drag coefficients for 
the aircraft. Results are printed out in tabular form and 
are graphed using a CALCOMP plotter. 

B. AIRCRAFT POWER REQUIREMENTS 

The heart of the program is the set of engine perfor- 
mance curves shown in Fig. 2. On the sea level side of the 
chart the curves are a linearized version of the original 
curves (Continental Motors Corporation, 195$) , and are a 




Figure 1 

General Program Flow Chart 
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good approximation to the originals. In the recommended 
cruising range the maximum error in the linearized curves 
is approximately 1 to 1.5$ in the mid-range RPM's and mani- 
fold pressures (or a maximum of approximately 2.5 brake 
horsepower), and much lower in the upper and lower ranges, 
where the original curves are very nearly linear. The 
altitude curves are a copy of the original curves. 

The curves are utilized in the following manner. For 
each RPM line in the sea level set of curves, the slope of 
the line, maximum manifold pressure, and maximum brake horse- 
power have been determined. For the altitude set of curves 
each RPM line has had slope and maximum brake horsepower 
determined. A question may arise at this point as to why, 
for a given RPM, the maximum brake horsepower is not the 
same for both curves and why, therefore, one horsepower 
value is not sufficient for any calculations needed. The 
answer is two-fold. First, in linearizing the sea level 
set of curves, the maximum points for the curves changed, 
whereas the maximum points for the altitude side did not; 
hence, two values. Secondly, since points on the RPM lines 
are determined using slopes and maximum points, it is 
necessary to keep both values so that calculations for both 
sets of curves will be accurate. 

A simplifying assumption was made concerning the mani- 
fold pressure lines on the altitude side. For purposes of 
this program the lines were assumed to be vertical. This 
introduces an error of not more than 1.1$, which yields at 
most a> difference of two brake horsepower. This assumption 
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Figure 2 

ENGINE PERFORMANCE CURVES 
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greatly simplifies the method by which manifold pressure is 
accounted for. Since the altitude curves are graphed as 
a linear function of the density ratio, a simple conver- 
sion is made from manifold pressure to density ratio and 
this converted ratio (called MAPSIG) is used for calcula- 
tions, as will be shown later. Table I shows the data used 
for the MAPSIG conversion. This is a purely mathematical 
manipulation, used for simplifying the equations by restating 
one variable in terms of another. This converted ratio is 
related to the critical operating altitude and is there- 
fore differentiated from the true density ratio (SIGMA) for 
the actual operating altitude, which is determined in the 
program from the flight altitude. Table II shows the data 
used for determining the flight altitude density ratio. 

Once the engine curve parameters have been determined, power 
requirements are calculated for a single engine. 

On the sea level side, brake horsepower (SSLBHP) is 
calculated as a function of the maximum brake horsepower 
(MAXBHP) and maximum manifold pressure (MAXMAP) for the 
operating RPM, slope of the RPM line (RSLOPE) , and operating 
manifold pressure (MAP) according to the equation: 

SSLBHP = MAXBHP - RSLOPE (MAXMAP - MAP) (l) 

The procedure for determining a value for brake horsepower 
at altitude is reasonably straightforward. The critical 
altitude (BHPCR) for a given manifold pressure is that 
altitude represented by the intersection of the operating 
RPM and manifold pressure lines. The converted density 
ratio (MAPSIG) discussed earlier defines this altitude. 
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TABLE I 



MAP SI G CONVERSION DATA 



MAP 


MAPSIG 


MAP 


MAPSIG 


MAP 


MAPSIG 


14. 


.545 


20. 


.730 


26 . 


.905 


15. 


.5 So 


21. 


.760 


27. 


.935 


l6. 


.615 


22. 


.790 


23. 


.965 


17. 


.640 


23. 


.320 


29. 


.990 


13. 


.675 


24. 


.355 






19. 


.705 


25. 


.330 







ALT 


TABLE II . 

ALTITUDE/DENSITY RATIO CONVERSION DATA 


SIGMA 


SIGMA 


ALT 


SIGMA 


ALT 


0 . 


1. 00000 


6000. 


,83586 


12000. 


.69317 


1000. 


.97106 


7000. 


.31064 


13000. 


.67133 


2000. 


.94277 


Sooo. 


.73601 


14000. 


.65002 


3000. 


.91512 


9000. 


.76196 


15000. 


.62923 


4000. 


,888 08 


10000. 


.73343 






5000. 


, 8 6167 


11000. 


.71555 
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A straight line is drawn connecting this point with the 
sea level brake horsepower value. Then a vertical line is 
drawn upward from the operating density ratio on the abscissa. 
If this line intersects the SSLBHP-BHPCR line, then that 
intersection defines the brake horsepower required for 
level flight at the operating altitude. If the vertical 
line lies to the right of the critically altitude then its 
intersection with the operating RPM line defines the 
required brake horsepower. Figure 3 shows an example of 
engine curve utilization. 

Critical brake horsepower (BHPCR) is calculated as a 
function of maximum brake horsepower for the operating RPM 
(MABHP) , the slope of the operating RPM line (ARSLOP) and 
the converted density ratio (MAPSIG) according to the 
equation: 

BHPCR = MABHP - ARSLOP (1.0 - MAPSIG) (2) 

The slope of the line connecting sea level and critical 
brake horsepowers (LSLOPE) is also a function of the con- 
verted density ratio: 

LSLOPE = (BHPCR - SSLBHP) / (1.0 - MAPSIG) (3) 

The operating altitude is compared with the critical altitude. 
If the operating altitude is the lower one, brake horsepower 
is calculated on the operating RPM line as: 

BHPALT = SSLBHP + LSLOPE (1.0 - SIGMA) (A) 
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If the operating altitude is above the critical altitude 
then brake horsepower is calculated on the operating RPM 
line as: 



BHPALT = MABHP - ARSLOP (1.0 - SIGMA) (5) 

It should be noted at this point that the actual flight 
data are based on thrust horsepower rather than brake horse- 
power. In the classroom the student takes an ideal brake 
horsepower, applies the necessary corrections and then mul- 
tiplies that value by propeller efficiency to get the actual 
thrust horsepower available. This program does exactly the 
opposite. From in-flight RPM and manifold pressure it cal- 
culates the actual thrust horsepower required. This is 
corrected, and then converted to shaft horsepower (it is 
assumed that shaft and brake horsepower are equal) by divid- 
ing by propeller efficiency. 

Since the engine performance curves are plotted in 
terms of brake horsepower, the nomenclature of the program 
up to this point has been written in terms of brake horse- 
power. However, since the horsepower value just found using 
equation (5) is actually a thrust horsepower, from here on 
it will be referred to as a thrust horsepower, until it is 
converted to shaft horsepower. 

Once the basic thrust horsepower value has been deter- 
mined, several correction factors are applied, the first 
of which is for temperature variation. The engine curves 
are plotted for standard temperature; a correction factor 
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ENGINE CURVE UTILIZATION 
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is supplied with them. First the temperature is changed 
from degrees Fahrenheit (TEMPF) to degrees Centigrade 
(DEGC) and the difference between operating and standard 
temperatures (TDIF ) is determined: 

DEGC = 5.0 (TEMPF - 32.0) / 9.0 (6) 

TDIF = 15.0 - DEGC \ (7) 

Then, using the correction factor, the single engine thrust 
horsepower corrected for temperature variation (CBHP) is: 

CBHP = BHPALT (1.0 -0.0016667 (TDIF)) (g) 

The total thrust horsepower (TBHPA) is found by 
multiplying the temperature corrected value by the number 
of engines: 

TBHPA = CBHP (NOENG) (9) 

Total shaft horsepower required at altitude is found 
by dividing thrust horsepower by propeller efficiency: 

TSHPA = TBHPA / PE (10) 

Propeller efficiency has been arbitrarily chosen as 0.75 
and is entered in the program using a DATA statement for 
ease of access, should subsequent evaluation dictate a 
change in the value. 

The total shaft horsepower is good for one altitude 
and one temperature. To make this information more useful 
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it is converted to an equivalent shaft horsepower (SSLPWR) 
which may be used to find horsepower at any altitude for 
the given temperature: 



SSLPWR = TSHPA /SIGMA (11) 

The relationship between velocity and power and their equiv 
alent counterparts is most easily shown by the following 
sequence of equations: 

1 -jj-2 

^ssl “ 2 

and 

1 tt2 

lalt = » P v 

also 

O' = /»//% > = O' 

SO 

lalt = * o' 1,2 

Define equivalent velocity: V e = V 

then 

1 TT 2 

^alt ~ 2 /^o e 

P = TV = DV 

= * S C d (O v 2 v = * s C d ,0 v 3 
= * S C d CT V 3 

Define equivalent power: P g = P ycr 

p /o=” = * S c d/ o o ai v 3 = * s c d /^ ( ycT V) 3 

= 4- s c, v 3 

d 'o .e 
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Once equivalent power has been calculated, the remaining 
step in the power requirements section is to change velocity 
to equivalent velocity (EQVEL): 

EQVEL = VEL J SIGMA (12) 

C. LIFT AND DRAG 

Lift is a function of density, velocity, characteristic 
area and the lift coefficient. In order to avoid the 
necessity of inputing or calculating a density for the 
flight altitude, the equations for determining lift and 
drag coefficients are based on equivalent velocity, reducing 
density to a constant (sea level value). 

Lift = weight = ■§• /& V 2 S 

= i /% or V 2 S C 1 

= i {Jv V ) 2 s c 1 
= £ /& ? 2 s c, 

* s ' o e 1 

For the calculations, equivalent velocity is changed from 
knots to feet per second (VEFPS): 

VEFPS = EQVEL (6000.) / 3600. (13) 

Then the lift coefficient (CLIFT) is found as a function 
of weight ( WT ) and equivalent velocity: 

CLIFT = 2 (WT) / (.0023769) (S) (VEFPS) 2 (14) 

The characteristic area, S, is defined by a data state- 
ment at the beginning of the program for ease of access. 

There are two values given for S (NASA TN-D-623&, 1971) o 
In this program 200 square feet is used, which includes tip 
tanks. Without them the wing area is 179 square feet. 
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Drag is a function of density, velocity, characteristic 
area and the drag coefficient. It is also a function of 
velocity and power: 

Drag = P / V = P e / V e 

= ^ V 2 SC d 

or equivalently, P e = ■§• /O S V e 

= i Po s ° d v e 3 

For the calculations, equivalent power is changed from 
shaft horsepower to foot-pounds per second (PEFPPS): 

PEFPPS = SSLPWR (550.) (15) 

Then the drag coefficient (CDRAG) is found as a function 
of equivalent power and equivalent velocity: 

CDRAG = 2 (PEFPPS) / (.0023769) (VEFPS) 3 (l6) 
Once all computations have been completed, a table of 
velocity, manifold pressure, brake horsepower at altitude, 
equivalent velocity, equivalent shaft horsepower, lift 
coefficient and drag coefficient is printed out. Then two 
calls are made to DRAW, one of IBM's plotting subroutines 
used in conjunction with a CALCOMP plotter. The first call 
is for a graph of equivalent power versus equivalent velo- 
city; the second is for a graph of lift coefficient versus 
drag coefficient. 

D. INTERPOLATION SUBROUTINE 

As stated earlier, the heart of the program is the 
data obtained from engine curves for the Cessna. Originally 
it was decided that these data would be used to formulate 
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equations with RPM as the independent variable. Then, 
given an operating RPM, all necessary information could be 
obtained using the equations, eliminating the need for 
table searches and interpolation. Several subroutines 
from IBM's Scientific Subroutine Package, including LSQPL2, 
CHBFT, and SPLIN, were utilized in an effort to formulate 
equations which would accurately represent the engine data. 
However, due to problems with convergence, extended pre- 
cision requirements for curve coefficients, and decreasing 
accuracy with increasing curve order, all of these routines 
proved unsatisfactory in producing the required equations. 

It was therefore decided to write a simple subroutine which 
would interpolate in a table of data from the engine curves. 

Subroutine LINT performs linear interpolation or 
extrapolation. Inputs to the subroutine are two columns 
of data and an independent variable for which LINT will 
interpolate the dependent variable. The data column 
corresponding to the independent variable must be in 
monotonically increasing order but may be arbitrarily 
spaced. Since the engine curves are for a single engine, 
flight with two engines allows operation of each engine at 
values below those shown on the engine performance curves, 
but not above. Subroutine LINT was therefore written to 
extrapolate below the limits of the curves but not above. 

In the event that the independent variable falls above 
the range of data given, an error message so stating this 
is printed out. 
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The method used by LINT is very simple. First a check 
is made to see if the variable is below the range for inter- 
polation; if this is so, an extrapolation is performed. 
Otherwise the column is searched until bracketing values 
are found, at which time an interpolation is performed. If 
a bracketing pair cannot be found, the variable is above 
the data range and the error message is printed. 

E. SORTING SUBROUTINE 

To make the program easier to use, it was written so 
that the student's data may be loaded in any order. It is 
not necessary, for example, to put the data in order of 
increasing velocity. Subroutine SORT does this ordering 
during program execution so that the printed output data 
is ordered in terms of increasing velocity. SORT is 
written to evaluate one column of data, and to put two 
associated columns of data in the same chronological order 
as the one being evaluated. 

The method used by SORT is simple. A logical flag 
is set equal to FALSE. Then sequential pairs of values 
in the first data column are compared. If the second 
value in any one of these pairs is less than the first 
value in the pair, the positions of the two values, and 
of the associated values in the other two columns, are 
switched and the logical flag is set equal to TRUE. 

When the entire column has been searched, the flag is 
checked; if it has been reset (to TRUE), this indicates 
the possibility that there may be more switches necessary 
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so the flag is reset to FALSE and the procedure is repeated. 
When a complete column search is made without resetting 
the flag, the data are in proper order and control returns 
to the main program. Like the main program, the subroutine 
is limited to 100 data points per column. 

III. PROGRAM USAGg, 



A. DATA LOADING 

The first data to be loaded following execution of 
the program are the engine data, consisting of the slope, 
maximum manifold pressure and maximum brake horsepower 
(again, brake horsepower values are used here because the 
engine data comes from curves which are plotted in terms 
of brake horsepower) for sea level performance, and slope 
and maximum brake horsepower for altitude performance, for 
each of nine sets of RPM lines. These data are shown in 
Table III. Following the permanent data are the student's 
in-flight data. 

The order in which in-flight data is loaded varies. 

The program gives the option of constant or variable RPM, 
and the input order depends on which option is selected. 

To simplify the student data loading, input data formats 
for both options are included in the program. These formats 
are similar to those used by James Melsa in his book of 
computer programs for linear control theory. They give 
card by card instructions for loading the data, including 
card number, column numbers, data description, and input 
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TABLE III 



ENGINE PERFORMANCE DATA 



RPM 


RSLOPE 
( BHP/lNCH) 


MAXMAP 
(INCHES HG) 


MAXBHP 

(HP) 


ARSLOP 
( BHP/lNCH ) 


MABHP 

(HP) 


1300 


6.45 


29.75 


150.0 


173.91 


150.0 


1900 


7.16 


29.70 


164.5 


136.67 


163.0 


2000 


7.37 


29.65 


175.0 


200.00 


175.0 


2100 


8.17 


29.60 


193.0 


216.67 


191.0 


2200 


£.72 


29.55 


203.5 


235.16 


207.5 


2300 


9.24 


29.50 


223.5 


254.94 


223.5 


2400 


10.00 


29.40 


239.0 


267.03 


237.5 


2500 


10 .56 


29.40 


243.0 


231.32 


249.0 


2625 


11.09 


29.35 


260.0 


294.50 


260.0 



RPM 


MAP 

(INCHES HG) 


TABLE IV 

EXAMPLE PROBLEM 

VEL RPM 

(KNOTS) 


DATA 

MAP 

(INCHES HG) 


VEL 

(KNOTS) 


2100 


17.5 


35.0 


2450 


20.0 


151.5 


2100 


19.0 


100.0 


2300 


13.0 


135.0 


2100 


19.75 


115.0 


2100 


13.0 


119.5 


2100 


20.3 


130.0 


2100 


13.0 


103.0 


2200 


22.0 


145.0 


2100 


15.0 


37.5 


2300 


23.6 


160.0 


2100 


16.0 


30.0 


2300 


23.6 


160.0 


2100 


18.0 


39.0 


2200 


22.0 


145.0 


2100 


13.5 


111.5 


2100 


20.1 


130.0 


2100 


19.0 


122.0 


2100 


19.2 


115.0 


2100 


21.0 


130.5 


2100 


18.6 


100.0 


2300 


21.5 


144.5 


2100 


17.5 


35.0 


2300 


23.0 


152.5 


2450 


23.0 


164.0 


2300 


24.0 


153.5 
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format. These formats are presented in the program in 
Appendix A. 

Following the input of constant flight parameters, 
such as altitude and weight, velocities and manifold pres- 
sures are read in. If RPM is variable, its values are 
read in; otherwise the constant RPM value is read in and 
duplicated for each set of velocity and pressure. This 
duplication is done for use in the sorting subroutine, 
which requires a column of RPM data, and for the printed 
output, which includes a column for RPM. Once all the 
data are read in and sorted, the engine curve performance 
parameters are interpolated for each RPM, if RPM is vari- 
able; if it is constant the parameters are evaluated only 
once, saving unnecessary computer time. Following this 
evaluation, the power, lift and drag requirements are 
calculated. 

B. GRAPHICAL OUTPUT 

The power requirements and lift and drag coefficients 
for the aircraft are graphed on the CALCOMP plotter, using 
IBM's DRAW subroutine. Each graph is eight inches by 
ten inches, with a title below. The abscissa and ordinate 
are automatically scaled by the subroutine and the scale 
information is printed below the graph. Data points are 
plotted on the graph as triangles. This was done for 
several reasons. If the subroutine is instructed to plot 
a continuous curve, it merely connects each point to the 
next with a straight line. This will produce a very 
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erratic curve in the case of data with any degree of spread, 
and v/here there are only a few data points to plot. The 
use of triangles allows the student to draw a smooth curve 
through the data. It also allows the data to be loaded in 
any order, rather than forcing the student, for example, to 
arrange his data in increasing order. This last point is 
admittedly academic, since addition of '•the sorting sub- 
routine eliminates any plotting problems which might have 
occurred from unordered data. 

IV. EXAMPLE PROBLEM 

Two actual sets of student data were combined for the 
example problem. RPM was variable. Table IV shows the RPM 
values and associated manifold pressures and velocities. 
Temperature varied slightly during the test flights; an 
average value was calculated for use in the program. Air- 
craft and flight parameters are shown in Table V. The 
order of the data deck is shown in Appendix B. 

Results of the program compared very favorably with 
hand calculated data. Table VI shows the tabulated output 
from the program. There was some spread in the data, which 
is easily visible in the sample graphical output shown in 
Appendix C . 
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TABLE V 



EXAMPLE PROBLEM - AIRCRAFT FLIGHT PARAMETERS 

Number of engines operating: 2 

Flight altitude (feet): 5000 

Temperature (degrees F): 76.55 

Aircraft weight (pounds): 4500 

Number of data points: 26 
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SAMPLE PROBLEM TABULATED RESULTS 
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VEL (I) RPM(I) MAP ( I ) BIIPALT(I) EQVEL(l) SSLPWR(l) CLIFT(l) CDRAG(l) 



V. PROJECT EVALUATION 



It is felt that the program has met the objectives 
previously set forth. A simple, easy-to-use program has 
been written which evaluates aircraft in-flight data and 
produces graphs of power versus velocil^y, and lift coeffi- 
cient versus drag coefficient. It yields results which are 
comparable in accuracy to results calculated by hand. 
Additionally, there is a considerable savings in man-hours 
required by the student to evaluate his data, which was the 
original impetus for this thesis. 

VI . AREAS FOR PROGRAM EXPANSION 

There are several possibilities for future expansion 
of the program. Temperature variation, though small, may 
have a noticeable effect on performance figures, and a pro- 
gram option allowing variable temperatures would be bene- 
ficial. The same is true for an option which would allow 
altitude to vary. Test flights are subject to weather, and 
clouds, for example, do not always allow a flight to be 
performed at a single altitude. Another variable factor is 
weight. Fuel accounts for nearly 16 percent of the aircraft 
maximum gross weight, and on a long flight this could 
appreciably alter performance data. An option which uti- 
lizes fuel readings might well be a worthwhile addition. 
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Another section which might be added to the program is 
one which would plot rate of climb curves for the aircraft. 
This could be a simple subroutine which would require inputs 
of time to climb through a given altitude change over a 
range of velocities. An easy conversion to feet per minute 
could be made, and the results plotted. 

Other future possibilities include the area of dynamic 
stability and control information, such as roll rates, which 
are covered in the second half of the flight evaluation 
techniques course sequence. 
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APPENDIX A 
COMPUTER PROGRAM 
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APPENDIX D 



PROGRAM GLOSSARY 
ALT - (input) Altitude (feet) 

ARSLOP - Slope of any altitude performance RPM line 

BHPALT - Single engine uncorrected thrust horsepower at 
altitude 

BHPCR - Critical horsepower at altitude for any given RPM 
and MAP 

Cl - Conversion factor for changing velocity to feet per 
second 

C2 - Constant for use in determining drag coefficient 

C3 - Constant for use in determining lift coefficient 

CBHP - BHPALT corrected for temperature variation 

CDRAG - Drag coefficient 

CLIFT - Lift coefficient 

DEGC - Temperature, degrees Centigrade 

DUMMY - Dummy variable used for temporary storage 

EQVEL - Sea level equivalent velocity (knots) 

IALT - Integer version of ALT 
IPTS — Integer version of POINTS 

ITB - Vector of constants for use in DRAW subroutine 

IWT - Integer version of VJT 

LINT - Linear interpolation subroutine 

LSLOPE - Slope of the line connecting sea level and altitude 
values of BHP ( SSLBHP and BHPCR) for any given RPM 
and MAP 

MABHP - Maximum horsepower for any given altitude performance 
RPM line 

MAP - (input) Manifold pressure (inches Hg) 
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MAPSIG - Performance curve equivalent value of density ratio 
for any given value of MAP 

MAXBHP - Maximum horsepower for any given sea level perfor- 
mance RPM line 

MAXMAP - Maximum manifold pressure for any given sea level 
performance RPM line 

NOENG - (input) Number of engines operating 
ONERPM - (input) Const ant/variable RPM flag 
PE - Propeller efficiency 

PEFPPS - Sea level equivalent power (foot pounds per second) 
POINTS - (input) Number of data points 

RA - Vector of altitudes used for converting altitude to 
density ratio 

RC - Matrix of engine performance curve data 

RC1-RC6 - Vector equivalents of data contained in RC 

RMAP - Vector of manifold pressures used for MAPSIG con- 
version 

RPM - (input) Engine revolutions per minute 

RS - Vector of density ratios used for converting altitude 
to density ratio 

RSIG - Vector of density ratios used for MAPSIG conversion 
RSLOPE - Slope of any sea level performance RPM line 
RTB1 - Vector of data used in DRAW subroutine 
RTB2 - Vector of data used in DRAW subroutine 
S - Aircraft characteristic area (square feet) 

SIGMA - Density ratio for a given altitude 
SORT - Sorting subroutine 

SSLBHP - Sea level power for any given RPM and MAP 

SSLPW - Sea level equivalent total shaft horsepower 

TBHPA - Total horsepower at altitude corrected for tempera- 
ture and number of engines 
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TDIF - Difference between standard temperature and tem- 
perature at flight altitude 

TEMPF - (input) Temperature at altitude (degrees Fahrenheit) 
TITLE1 - Graph title information for DRAV/P subroutine 
TITLE 2 - Graph title information for DRAW subroutine 
TSHPA - TBHPA corrected for propeller efficiency 
VEFPS - Sea level equivalent velocity (feet per second) 

VEL - (input) Velocity (knots) 

VJT - (input) Aircraft weight (pounds) 
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